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ABSTRACT 

We discuss the lifetimes and evolution of clumps and cores formed as turbulent density fluctuations 
in nearly isothermal molecular clouds. In order to maintain a broad perspective, we consider both the 
magnetic and non-magnetic cases. In the latter, we argue that clumps are unlikely to reach a hydrostatic 
state if molecular clouds can in general be described as single-phase media with an effective polytropic 
exponent 7 C < 4/3. In this case, clumps are expected to be short-lived, either proceeding directly to 
collapse, or else "rebounding" towards the mean pressure and density as the parent cloud. Rebounding 
clumps are delayed in their re-expansion by their self-gravity. From a simple Virial Theorem calculation, 
we find re-expansion times ~ 1 Myr, i.e., of the order of a few local free-fall times tf c = Lj iC /c, where 
Lj )C = y/ irc 2 /Gp c is the Jeans length at the clump's mean density p c and c is the isothermal sound 
speed. 

In the magnetic case, we present a series of driven-turbulence, ideal-MHD numerical simulations in 
which we follow the evolution of clumps and cores in relation to the magnetic criticality of their "parent 
clouds" (the numerical boxes). In subcritical boxes, magnetostatic clumps do not form. A minority of 
moderately-gravitationally bound clumps form which however are dispersed by the turbulence in < 1.3 
Myr. An estimate of the ambipolar diffusion (AD) time scale at the physical conditions of these cores 
gives characteristic times > 1.3 Myr, suggesting that these few longer-lived cores can marginally be 
"captured" by AD to increase their mass-to-flux ratio and eventually collapse, although on time scales 
not significantly longer than the dynamical ones. In supercritical boxes, some cores manage to become 
locally supercritical and collapse in typical time scales of 2 tf c (~ 1 Myr). In the most supercritical 
simulation, a few longer-lived cores are observed, which last for up to ~ 3 Myr, but these end up re- 
expanding rather than collapsing, because they are sub- Jeans in spite of being super-critical. Fewer 
clumps and cores form in these simulations than in their non-magnetic counterpart. 

Our results suggest that a) Not all cores observed in molecular clouds will necessarily form stars, and 
that a class of "failed cores" should exist, which will eventually re-disperse, and which may be related to 
the observed starless cores, b) Cores may be out-of-equilibrium, transient structures, rather than quasi- 
magnetostatic configurations, c) The magnetic field may help reduce the star formation efficiency by 
reducing the probability of core formation, rather than by significantly delaying the collapse of individual 
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1. INTRODUCTION 

One of the most important goals in the study of star for- 
mation is to understand the state and physical conditions 
of the molecular cloud cores from which the stars form. 
The prevailing view (which we hereafter refer to as the 
"standard (magnetic support) model" of star formation; 
e.g., Mouschovias 1976a, b; Shu, Adams & Lizano 1987) is 
that low-mass star-forming clumps evolve into dense cores 
along a sequence of quasi- magnetostatic states as magnetic 
support is lost by ambipolar diffusion. This process has 
a relatively long characteristic time scale, typically esti- 
mated to be significantly larger than the local clump free- 
fall time (e.g., Shu et al. 1987; McKee et al. 1993). Dur- 
ing this time, the contracting clump is supported against 
its self-gravity by the magnetic field in the direction per- 



pendicular to it, and by a combination of thermal and 
micro-turbulent pressures along the field (e.g., Lizano & 
Shu 1989; McLaughlin & Pudritz 1996). 

The support from magnetic fields is included through 
the consideration of the mass-to-magnetic flux ratio of 
the core. For isolated or periodic structures, in which 
mass cannot be accreted from the surroundings, and un- 
der ideal MHD conditions, the magnetic flux is conserved, 
and therefore so is the mass-to-flux ratio (Chandrasekhar 
& Fermi 1953; Mestel & Spitzer 1956). Thus, self-gravity 
cannot overcome the magnetic support if the mass-to-flux 
ratio is smaller than some critical value, and collapse can 
only occur as the neutrals in the gas slowly slip across field 
lines by ambipolar diffusion (see, e.g., Mestel & Spitzer 
1956; Mouschovias & Spitzer 1976; see also the reviews by 
Shu et al. 1987 and McKee et al. 1993). 
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Additionally, although in the standard model the mag- 
netic support is fundamental, discussions of non-magnetic 
equilibrium configurations are very often encountered in 
the literature. These configurations have generally ci- 
ther finite or diverging central densities. The former 
are generally based on the so-called Bonnor-Ebert (BE; 
Ebcrt 1955; Bonnor 1956) spheres, which are solutions of 
the Lane-Emden equation, and can be either truncated 
at some finite radius, beyond which a hotter, less dense 
medium is assumed to confine the sphere, or infinite ( "ex- 
tended"). Truncated configurations have recently been 
taken as "templates" for comparison with observations of 
column density profiles in real cores (e.g., Johnstone et al. 
2000; Alves, Lada & Lada 2001; Evans et al. 2001; Harvey 
et al. 2001), even though Shu et al. (1987, sec. 4.1) had 
already pointed out that such truncated structures may 
have little relevance to actual cores in molecular clouds be- 
cause they are in contact with a large reservoir of cold gas, 
and that a maximum contrast of 14 contradicts the obser- 
vational evidence that molecular clouds typically have H2 
densities ranging from 10 2 cm -3 (average) to well in excess 
of 10 6 cm~ 3 in cores (see, e.g. Goldsmith 1987.) Various 
geometries other than spherical have also been considered 
for non-magnetic equilibrium solutions (e.g., Curry 2000; 
Lombardi & Bertin 2001). Concerning structures with di- 
verging central densities, the best known example is the 
singular isothermal sphere (e.g., Shu 1977). 

On the other hand, it is well established that the molec- 
ular clouds within which the cores form arc turbulent, with 
lincwidths that are supersonic for scales > 0.05 pc (e.g., 
Zuckerman & Evans 1974; Larson 1981; see also the review 
by Blitz 1993 and references therein), generally interpreted 
as magnetohydrodynamic (MHD) supersonic turbulent 
motions. Being supersonic (and possibly super-Alfvenic; 
see below), these motions in general dominate the sup- 
port against the clouds' self-gravity, with thermal pres- 
sure being dominant only at subsonic scales (Padoan 1995; 
Vazquez-Semadeni, Ballesteros-Paredes & Klesscn 2003a). 
In this environment, the molecular cloud clumps and cores 
are in fact likely to be the turbulent density fluctuations 
themselves within the clouds (von Weizsacker 1951; Bania 
& Lyon 1980; Scalo 1987; Elmegreen 1993; Ballesteros- 
Paredes, Vazquez-Semadeni & Scalo 1999;Klessen, Heitsch 
& Mac Low 2000; Padoan et al. 2001; Heitsch, Mac Low 
& Klessen 2001). 

Moreover, in a turbulent environment, the density and 
magnetic field strength are only very poorly correlated lo- 
cally, at least when gravity is not important in the forma- 
tion of the density sturctures (Passot, Vazquez-Semadeni 
& Pouquet 1995; Padoan & Nordlund 1999; Ostriker, 
Stone & Gammie 2001; Kim, Balsara & Mac Low 2001; 
Passot & Vazquez-Semadeni 2003). Thus, the clumps and 
cores can in principle be formed with a distribution of 
mass-to-flux ratios. For ideal isolated or periodic config- 
urations, this distribution is bounded from above by the 
mass-to-flux ratio of the parent cloud (cf. §3.3) for as 
long as ambipolar diffusion is negligible (i.e., at densities 
< 10 4 cm' 3 ). Real molecular clouds are not isolated nor 
periodic, and can accrete mass along the magnetic field to 
eventually become supercritical (Hartmann, Ballesteros- 
Paredes & Bergin 2001), although this is likely to occur 
on time scales much longer than those of core evolution. 



So, for the purpose of studying core evolution, the mass-to- 
flux ratio of the parent cloud can be considered constant. 

Within this framework, if the parent cloud is subcritical, 
then in principle supercritical cores can only form through 
the action of ambipolar diffusion, in agreement with the 
scenario of the standard model. Indeed, numerical simula- 
tions of supersonic turbulence in subcritical boxes without 
ambipolar diffusion have shown that collapse does not oc- 
cur (Ostriker, Gammie & Stone 1999; Heitsch et al. 2001). 
On the other hand, if the parent cloud is supercritical, 
then in principle both sub- and supercritical substructures 
can form, even under ideal MHD. The latter can collapse 
if they are also gravitationally unstable with respect to 
the thermal pressure, a case we will sometimes refer to as 
"super-Jeans" . 

Observationally, several recent studies (e.g., Crutcher 
1999, 2004; Bourke et al. 2001; Crutcher, Heiles & Troland 
2002) suggest that molecular clouds are the objects with 
highest mass-to-flux ratios (within the uncertainties) in 
the hierarchy going from diffuse medium to molecular 
cloud cores. Furthermore, Nakano (1998) has suggested 
that the large observed column densities of cores are dif- 
ficult to reconcile with the relatively low column density 
enhancements expected for subcritical cores with respect 
to their parent clouds, while Padoan and coworkers (e.g., 
Padoan & Nordlund 1999; Padoan at al. 1999; Padoan, 
Goodman & Juvela 2003) have performed several synthetic 
observations of cores in simulations of MHD turbulence 
suggesting that molecular clouds are super-Alfvenic. Fi- 
nally, like their parent molecular clouds, the clumps and 
cores (density peaks) within them are not isolated, but are 
connected to their surrondings, and so their masses are 
not fixed, especially in a turbulent environment. The non- 
isolated nature of the cores has been taken into account 
within the standard model (e.g., Lizano & Shu 1989; Curry 
2000), but still under the assumption of a quasi-static evo- 
lution. 

There thus appears to exist a conceptual gap between 
the standard assumption of quasi- hydrostatic clump/core 
evolution and the turbulent nature of the clouds in which 
they form, the natural question being whether cores 
formed dynamically in a turbulent medium can settle into 
the (magneto-)hydrostatic equilibria assumed in the stan- 
dard model. Instead, they might as well be either pushed 
directly into collapse, or re-expand and merge back with 
their environment. Recently, Li & Nakamura (2004) have 
addressed this issue by considering the formation of dense 
cores in two-dimensional simulations of sheet-like, strongly 
magnetized clouds including a prescription for ambipolar 
diffusion. In the present paper we consider the problem 
from a more general point of view, by discussing clump and 
core formation in both strongly and weakly (including non- 
magnetic) magnetized clouds, investigating the lifetimes 
of clumps formed by transient turbulent fluctuations and 
the evolution of their Jeans number and magnetic criti- 
cality in relation to that of their parent molecular clouds. 
We begin (§2) by presenting and discussing the main as- 
sumptions that underlie our approach, as well as obser- 
vational and theoretical evidence that supports them. We 
then proceed to give some general considerations in (§3) in 
both the non-magnetic and magnetic cases. With respect 
to the former, in §3.1 we present a qualitative discussion 
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of the reasons why hydrostatic configurations are unlikely, 
in single-phase, nearly isothermal, non-magnetic turbulent 
flows. This implies that dynamically compressed regions 
must either proceed to collapse right away or re-expand. 
In §3.2 we then give a simple estimate of the re-expansion 
time for clumps that do not collapse (which we call "failed" 
clumps). In the magnetic case, in §3.3 we give a qualitative 
discussion of the formation and nature of cores in relation 
to the magnetic criticality of their parent cloud. 

In §§4 and 5 we then present a series of three- 
dimensional numerical simulations aimed at studying the 
lifetimes of cores in non-magnetic and magnetic sub- and 
supercritical environments, first describing the numerical 
method and parameters (§4.1), resolution considerations 
(§4.2) and core analysis procedure (§4.3), and then the re- 
sults (§5). Finally, in §6 we discuss some implications and 
caveats of our results, and in §7 we give a summary and 
concluding remarks. 

2. ASSUMPTIONS 

In this section, we describe the assumptions underlying 
the discussions and calculations in this paper, as well as 
the evidence supporting them. Nevertheless, as will be 
clear from the discussion, the assumptions are not conclu- 
sively proven, and therefore it must be kept in mind that 
our results are of course subject to the applicability of our 
assumptions. 

2.1. Single-phase isothermal medium 

Throughout this paper we assume that the gas within 
molecular clouds responds isothermally to compressions 
(i.e., that it obeys an isothermal "effective" equation of 
state P = c 2 p, where P is the thermal pressure, p is the 
density, and c is the sound speed), and that it is at the 
same temperature throughout the volume considered. Al- 
though this is the most common assumption for numeri- 
cal simulations of molecular clouds (see, e.g., the work of 
Padoan, E. Ostriker, Mac Low, Klessen, Bate, Bonnell and 
their coworkers; see also the reviews by Vazquez-Semadeni 
et al. 2000; Mac Low & Klessen 2004) because molecular 
gas is typically at temperatures ~ 10-15 K, the true ther- 
modynamic behavior of molecular gas is still an open issue. 

The cores in the survey of Jijina, Myers & Adams (1999) 
span almost an order of magnitude in temperature, from 
~ 5 to 40 K, although indeed with a strong concentration 
between 10-15 K, and presumably owing in most cases 
to direct heating by nearby stellar objects, rather than to 
PdV mechanical heating. In the absence of direct stel- 
lar /protostellar heating, the gas tends to have more uni- 
form temperatures. Indeed, the subsample of cores with no 
embedded clusters of Jijina et al. has a narrower tempera- 
ture distribution, while a study of the dark cloud TMC-1 
by Pratap et al. (1997) gives nearly uniform temperatures, 
between 8 and 10 K, even though the density spans nearly 
an order of magnitude. 

Nevertheless, precise isothermality is not expected in 
molecular clouds. Scalo et al. (1998) summarized results 
on the likely values of the exponent 7 e that appears in 
polytropic equations of state of the form P oc /? 7 ° (see also 
Myers 1978; Larson 1985). For molecular gas, Scalo et 
al. found 1/3 < 7c < 1 for 10 3 cm" 3 < n < 10 5 cm" 3 
in the absence of embedded protostellar sources. Recent 



calculations of the effective polytropic index by Spaans 
& Silk (2000) find a similar range for present-day Galac- 
tic molecular gas. More recently, Li, Klessen & Mac Low 
(2003) have investigated the dependence of the fragmenta- 
tion process fragmentation on the value of the polytropic 
exponent in molecular clouds. 

For the purposes of the present paper, the most im- 
portant thermodynamic consequences of assuming isother- 
mality (i.e., 7e = 1) are that the increase in thermal pres- 
sure during gravitational collapse cannot eventually halt it 
(something which would require j e > 4/3 for spherically- 
symmetric collapse; e.g., Hayashi 1967; Larson 1969; Low 
& Lynden-Bell 1976; Masunaga & Inutsuka 1999; see also 
the discussion in §3.1), and that the medium consists of 
a single thermodynamic phase. A two-phase, or "ther- 
mally bistable" , nature would require the existence of an 
intermediate thermally unstable density range with 7 e < 
between density ranges with 7 C > (e.g., Field, Goldsmith 
& Habing 1969). However, we see that the 7 c -range esti- 
mated by Scalo et al. (1998) and Spaans & Silk (2000) for 
the densities of interest, 1/3 < j c < 1, is completely con- 
tained within the interval < 7 < 4/3, and therefore also 
does not allow for neither of those two phenomena. Thus, 
in this sense, our isothermal assumption is consistent with 
those 7 -range estimates. 

However, we must point out that Gildcn (1984) has sug- 
gested that a thermally unstable range with j c < may 
exist in molecular gas, and, in consequence, allow for the 
possibility of thermal bi-stability. If this does indeed oc- 
cur, then our discussion of the impossibility of producing 
Bonnor-Ebcrt-type hydrostatic structures in §3.1, based 
on the non-availability of a hotter, more tenuous confin- 
ing medium, is invalidated. Moreover, Blitz & Williams 
(1999) have argued that the HI gas may pervade the vol- 
ume within clumps, by showing an anticorrelation between 
the CO and the HI emission. This could also provide 
a confining medium for BE- type structures. Ultimately, 
observations and detailed studies of the heating/cooling 
processes in the gas should determine the goodness of the 
isothermal assumtion for molecular gas. 

Other relatively minor consequences of using an isother- 
mal prescription are that sound waves cannot steepen into 
shocks, and that the production of vorticity through the 
baroclininc term is inhibited. The latter limitation extends 
to polytropic prescriptions in general (see, e.g., Vazquez- 
Semadeni et al. 1996). This suggests that forthcoming 
improved molecular cloud simulations should solve the in- 
ternal energy equation, with reasonably realistic prescrip- 
tions for the cooling and heating processes within them, 
as has been done for simulations of the global ISM (see, 
e.g., the review by Vazquez-Semadeni 2002). 

2.2. Driven turbulence regime 

Another assumption underlying our work in this paper 
is that the turbulence in molecular clouds is driven, rather 
than decaying, as we have used almost exclusively simula- 
tions of driven turbulence. Observationally, this is moti- 
vated by the fact that turbulence is ubiquitous in molec- 
ular clouds (e.g., Larson 1981; Blitz 1993; Heyer & Brunt 
2004), including starless ones such as the so-called Mad- 
dalena's cloud (Maddalena & Thaddeus 1985). This sug- 
gests a universal origin and maintenance mechanism for 
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the clouds' turbulence (Heyer & Brunt 2004). It has been 
suggested that the turbulence in the clouds may arise from 
the very accumulation mechanism that forms them out 
of the more diffuse atomic gas (Ballesteros-Paredes et al. 
1999a; Koyama & Inutsuka 2002; Vazquez-Semadeni et al. 
2003a,b; see also Hunter et al. 1986; Waldcr & Folini 1998; 
Vazquez-Semadeni 2003), through dynamical instabilities 
in the compressed regions. In this case, the turbulence 
in molecular clouds is driven for as long as the accumu- 
lation process continues, while they probably disperse af- 
terwards, as in the recent simulations by Clark & Bonnell 
(2004). 

In this scenario, the clouds are in fact part of the turbu- 
lent cascade from the largest scales down to the dissipative 
scales in the ISM, and the driving therefore occurs from the 
large scales. This is supported by recent comparisons be- 
tween observations and either numerical simulations (Os- 
senkopf & Mac Low 2002) or fractional-Brownian-motion 
velocity fields (Brunt 2003), in which the observed fields 
are most consistent with the large-scale-driven turbulence 
numerical cases. Moreover, using Principal Component 
Analysis, Brunt & Heyer (2002) have estimated the spec- 
tral index of the turbulence in 23 fields of the FCRAO CO 
survey of the Outer Galaxy, finding that it is most con- 
sistent with the index found in numerical simulations of 
driven turbulence. Thus, we assume that the appropriate 
model for molecular clouds is that of driven turbulence, 
with injection at the largest scales. 

In this paper we use the standard prescription of driv- 
ing the flow in Fourier space. Although this allows precise 
control of spectral properties of the turbulence, such as the 
driving scale, it has the downside of injecting the energy 
everywhere in physical space, and is therefore not the most 
realistic energy injection mechanism. More realistic sim- 
ulations of molecular clouds should include physical driv- 
ing mechanisms, which, according to our discussion here, 
amount possibly to the cloud formation process itself. In 
any case, the driving scheme we use here should not di- 
rectly interfere with the cores' evolution, as it is applied 
on scales at least 10 times larger than typical core sizes. 
Specifically, the driving is applied at scales of 1/2 the box 
size, while clumps and cores in the simulations have typ- 
ical sizes 0.08 and 0.04 times the box size (0.3 and 0.15 
pc), respectively. 

2.3. Neglect of ambipolar diffusion 

The numerical simulations we present in §5 are of ideal 
MHD turbulence (subject only to the intrinsic numerical 
diffusion of the numerical scheme), and do not include 
ambipolar diffusion (AD). This implies that, in the case 
of subcritical simulations, gravitational collapse is com- 
pletely suppresed. In this regard, the recent simulations 
by Li & Nakamura (2004) including AD, give a more ac- 
curate description of the effect of AD. On the other hand, 
those authors restricted their simulaions to two dimensions 
and did not consider the possibility of supercritical clouds, 

1 Note that for the realization of two- and one-dimensional equilibria (i.e., with cylindrical and plane-parallel geometries, respectively), the 
values of the critical polytropic exponent are 7 C = 1 and -y c = 0, respectively (Ostriker 1964, 1965; Larson 1985; Vazquez-Semadeni et al. 
1996). However, for the three-dimensional, non-magnetic case we are considering here, we do not expect such configurations to be applicable, 
as nothing locally restricts the flow to a lower dimensionality. For example, filament production is often observed in simulations, but collapse 
then proceeds along the filament (e.g., Ballesteros-Paredes 2004; Klessen & Ballesteros-Paredes in preparation), not perpendicularly to it as 
would correspond to cylindrical geometry. On the other hand, in the presence of a magnetic field, an anisotropy arises which tends to make 
the turbulent compressions more one-dimensional. This is why stable magnctostatic equilibria can exist in a sufficiently strongly magnetized 
case, even though contraction can proceed freely along field lines. 



which we do here. 

As a zeroth-order approximation to determining the ef- 
fect of AD in our subcritical simulation, in §5.2.1 we esti- 
mate the classical AD time scale tad under the conditions 
of the clumps that form there, which can then be com- 
pared to the clumps' lifetimes. If tad is much larger than 
the latter, then AD is not expected to have time to act 
during the cores' duration. If it is comparable or smaller, 
then it can "capture" a transient, subcritical core, and 
increase its mass-to-flux ratio, eventually causing its be- 
coming supercritical and collapse, although on time scales 
comparable to the dynamical lifetimes of the clumps we 
find in the simulations, thus not significantly delaying the 
collapse time. 



3. GENERAL CONSIDERATIONS 
3.1. The non-magnetic case 

In this section we suggest that the formation of hy- 
drostatic configurations is unlikely in single-phase, nearly 
isothermal (or, more generally, with < j c < 4/3), non- 
magnetic flows. This is because, for nearly-isothermal 
flows, the production of significant density fluctuations re- 
quires supersonic turbulent compressions. These fluctua- 
tions are in general transient, and so the density enhance- 
ments they produce ("clumps") must also be so, unless 
the latter happen to "land" near a stable (attracting, in 
the language of dynamical systems) hydrostatic equilib- 
rium between self-gravity and thermal pressure (not tur- 
bulent pressure, because it is locally too highly variable 
for establishing a static configuration). In this case, after 
the external turbulent compression subsides, the core may 
evolve towards this equilibrium rather than towards either 
re-expansion or gravitational collapse. 

It follows that hydrostatic cores can occur only if there 
exist stable equilibria. The best known such configura- 
tions in the non-magnetic, isothermal case are the well- 
known Bonnor-Ebert (BE) solutions of the Lane-Emdcn 
equation. However, the stability of these configurations 
relies on the fact that they are truncated at a finite radius 
from their center, being confined by a hotter, negligible- 
density medium. The boundary between the cold, dense 
medium and the hot, low-density one amounts to a transi- 
tion between two stable phases (with j c > 0) mediated by 
a (non-populated) thermally unstable density range with 
7c < 0. 

Of course, a hot confining medium is not necessary in 
general. For polytropic situations with 7 C > 7 C , where 
7 C = 4/3 (Chandrasekhar 1961, §117; Yabushita 1968; 
McKce et al. 1993, sec. VI. C; Vazquez-Semadeni, Pas- 
sot & Pouquet 1996), spherically-symmetric equilibria can 
be found even for a vanishing external pressure. 1 This is 
why stars can be formed as stable entities from highly 
anisotropic, dynamic, time-dependent accretion (Hart- 
mann et al. 2001). However, for bounded systems with 
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7 C < 4/3, the boundary pressures are indispensible in es- 
tablishing stable equilibria. 2 Such boundaries arc, how- 
ever, ruled out in single-phase media, in which no abrupt 
transition to a hotter, more tenuous medium at the same 
pressure is possible. 

We must examine, then, the possibility of extended (i.e., 
non-truncated) stable equilibrium configurations, in which 
the core merges smoothly into its environment. To our 
knowledge, this remains an open question. The best- 
known extended equilibrium configuration is the singular 
isothermal sphere, which, however, is known to be un- 
stable, and has been used in fact as the initial condition 
for gravitational collapse calculations (e.g., Shu 1977). A 
more realistic (though still spherical) configuration would 
be one with a finite central density. Such a configura- 
tion can be obtained by numerically integrating the Lane- 
Emden equation to arbitrarily large radii. However, this 
implies that it would be equivalent to a BE sphere of arbi- 
trarily large radius, and would thus be gravitationally un- 
stable 3 (wc thank D. Galli for suggesting this argument). 

Other possibilities exist, however. Curry (2000) has 
found a new class of periodic equilibrium structures in 
cylindrical geometry, which he has proposed as candidates 
for the formation of cores in filamentary clouds. However, 
he did not give definitive results concerning the stability of 
these structures, and in fact concluded that full numerical 
simulations may be necessary to resolve the issue. 

In this respect, numerical simulations of non-magnetic 
self-gravitating polytropic turbulence with j c < 1 sug- 
gest that hydrostatic configurations do not form: the tur- 
bulent density fluctuations either re-expand, or proceed 
directly to collapse. This was first shown by Vazquez- 
Semadeni et al. (1996) for two-dimensional turbulent flows 
with 7o < 0.3, 4 in which many generations of clumps were 
seen to appear and disappear in roughly their crossing 
time, until finally a strong enough compressive event oc- 
curred as to produce a gravitationally unstable core, which 
then proceeded to collapse in roughly its own free-fall time. 
This result has subsequently continued to be found in 3D 
simulations of non-magnetic, isothermal (e.g., Klcssen et 
al. 2000) as well as polytropic (Bate, Bonnell & Bromm 
2002; Li, Klcssen & Mac Low 2003) turbulence, in which 
the formation of hydrostatic objects is in general not ob- 
served (R. Klessen, private communication) for j c < 4/3. 
The turbulent density fluctuations either re-expand or pro- 
ceed directly to collapse in their own free-fall time. 

Finally, it is worth noting that Clarke & Pringle (1997) 
have pointed out that cores cool mainly through optically 
thick lines, but are heated by cosmic rays, and therefore 
may be dynamically unstable, as velocity gradients may 
enhance local cooling. 

3.2. Re-expansion time of "failed" compressions 

The arguments above suggest that a significant fraction 
of the clumps and cores in molecular clouds may be headed 
towards re-expansion rather than collapse. This possibility 
is often overlooked in the literature, which has generally 
focused on the gravitational collapse of cores to form stars, 

2 It is interesting to note that Hunter (1977) found a class of unstable BE spheres that do not collapse, but instead undergo large-amplitude 
oscillations. However, these still rely on the presence of an external confining medium at constant pressure and neglible density. 

3 Again, in practice, the typical core-to-cloud mean density ratio of ~ 30-100 is well into the unstable regime for BE-spheres. 

4 In this 2D case, 7 cr = 1. 



and considers starless cores as "pre-stellar" (with a few no- 
table exceptions; e.g., Taylor, Morata & Williams 1996). 
However, in the scenario advocated in this paper, a frac- 
tion of the starless cores may never actually form stars, 
and instead rc-cxpand and disperse back into their parent 
molecular cloud, being natural candidates to correspond 
to at least some of the observed starless cores. 

In the absence of self-gravity, one would expect that the 
re-expansion of a spherical core should be nearly identi- 
cal to the compression that produced it, run backwards 
in time. However, in the self-gravitating case, the re- 
expansion process must occur on a longer time scale than 
the compression because of the retarding action of self- 
gravity. It is thus of interest to estimate their expected 
lifetimes. 

A crude estimate of the re-expansion time can be given 
in terms of the Virial Theorem (VT) , similarly to the anal- 
ysis performed by Hunter (1977, sec. III.) to investigate 
the classes of solutions of unstable BE spheres. We ex- 
tend his treatment to estimate the re-expansion time scale, 
although we consider the case of zero external pressure, 
because a) the density contrast between cores and their 
parent molecular clouds are large enough (up to factors 
~ 100) that the external thermal pressure is negligible 
during the initial stages of the evolution with which we 
are concerned, and b) in our case of an extended core, 
with a smooth density profile, the whole profile varies in 
time, and therefore the pressure is also variable at any 
boundary that we choose to define (which, as a matter of 
fact, would be rather arbitrary). The full time-dependent 
problem is what is solved by the numerical simulations, 
but for the purpose of a simple estimate, here we consider 
the free re-expansion of a gas sphere subject exclusively to 
its self-gravity and internal pressure. 

The VT for an isothermal spherical gas mass ( "cloud" ) 
of volume V and mean density fi in empty space is 



\l = 3Mc 2 - aGM 2 /R, 



(1) 



where the overhead dots indicate time derivatives, M = 
J v polV is the cloud's mass, R is its radius, p is the den- 
sity, I = j v pr 2 dV is its moment of inertia, c is the sound 
speed, and a is a factor of order unity. We can obtain an 
evolution equation for the cloud's radius by replacing the 
radius-dependent density by its mean value in the expres- 
sion for I, to find I w MR 2 . Thus, 



±1, 
2 



M 



R(t)R(t)+R 2 (t) 
From equations (1) and (2), we obtain 



R(t)R(t) + R 2 (t) 



= 3c 2 



aGM/R. 



(2) 



(3) 



This equation can be integrated analytically, with solution 
V(r 2 - I) 2 - (ri - I) 2 + 



T = 



1 



In 



-l + V(r 2 -l) 2 -(r-i-l)' 



(4) 
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where ri and r 2 are the initial and final radii of expansion, 
normalized to the equilibrium radius R c — aGM/3c 2 , and 
t = t/tff is the time, non-dimensionalized to the free-fall 
time tff = Re/c. A characteristic re-expansion time can 
be defined as the time required to double the initial ra- 
dius (i.e., r 2 = 2n, a "2-folding" time), starting from an 
initial condition r\ > 1. Figure la shows this characteris- 
tic time as a function of n. We see that when n is very 
close to unity (i.e., linear perturbations from the equilib- 
rium radius), the re-expansion time can be up to several 
times the free-fall time. Values of n moderately larger 
than unity have the shortest re-expansion times, because 
the initial force imbalance is greater, yet the final size is 
still not much larger than twice the equilibrium radius. 
Finally, for larger initial radii (far from the equilibrium 
value), the re-expansion time approaches that of free ex- 
pansion at the sound speed. These times are also long, 
but not because of a low expansion velocity, but because 
the distance to be traveled by the clump boundary from n 
to r 2 increases proportionally to n. Figure lb shows an 
crude "observer's perspective" of the re-expansion time, 
giving the time the core's density stays above 10% of the 
equilibrium-configuration density p c . This emulates the 
time the core would be visible to a high-density tracer sen- 
sitive to densities > 0.1 p c . This perspective clearly shows 
that the longest re-expansion times of an object that can 
be called a "core" correspond to those that come closest 
to equilibrium. 

We conclude that the re-expansion time is at least larger 
than twice the free-fall time, making the probability of ob- 
serving a core in this process larger by this factor than that 
of observing a free-falling core. This is consistent with the 
fact that molecular clouds are generally observed to con- 
tain more starless than star- forming cores (e.g., Taylor et 
al. 1996; Lee & Myers 1999; see also Evans 1999 and ref- 
erences therein). 

A number of points are worth noting here. First, the 
database put together by Jijina et al. (1999), reports a 
larger fraction of star- forming than starless cores but, as 
those authors recognize in their section 7.1, their sample 
is not free of biases, with starless cores not being as well 
represented as the star-forming ones. This is due to the se- 
lection criteria used to select the cores in the surveys from 
which these authors put their database together. Instead, 
the Lee & Myers (1999) dataset was designed to be as free 
as possible from this bias, their cores being selected op- 
tically from the STScI Digital Sky Survey. Second, some 
observations, such as those of ammonia cores, may be bi- 
ased towards selecting collapsing cores, if the densities they 
probe are very large; i.e., if they probe densities not regu- 
larly reached by non-collapsing clumps. The cores in the 
Jinina et al. sample, which is an ammonia-core database, 
have larger typical number densities (n > 10 4 cm~ 3 ) and 
a smaller starless fraction than those in the Lee & Myers 
optically-selected sample (n~7x 10 3 cm -3 ). Indeed, in 
the simulations reported in §5, the failed clumps rarely 
exceed peak densities ~5x 10 4 cm~ 3 . 

3.3. The magnetic case 

In the magnetic case, the classical Virial-Theorem 
(VT) analysis (Chandrasekhar & Fermi 1953; Spitzer 
1968; Mouschovias 1976a, b; Mouschovias & Spitzer 1976; 



Zweibel 1990) predicts the existence of sub- and super- 
critical configurations (e.g., Shu et al. 1987) depending 
on whether the mass-to-magnetic flux ratio is below or 
above a critical value (M/<ft) c . Subcritical configurations 
are known not to be able to collapse gravitationally unless 
the magnetic flux is lost by some process such as ambipolar 
diffusion. Supercritical configurations, on the other hand, 
are analogous to the non-magnetic case, except for the fact 
that the cloud behaves as if having an "effective" mass, re- 
duced by an amount equal to the critical mass (which de- 
pends on the magnetic field strength) (e.g., Spitzer 1968, 
§11.3.b). 

A magnetized clump or core (we generically use the 
term "core" in this section) may thus be characterized 
by two nondimcnsional numbers, the first one being its 
Jeans number (J c = l/Lj, c ), where I is the characteristic 
core size and Lj, c is the Jeans length at the mean core 
density, measuring whether the core is gravitationally un- 
stable with respect to the thermal support (J c > 1, or 
"super- Jeans" ) or not. The other is the core's mass-to- 
magnctic-flux ratio expressed in units of the critical value 
for collapse, p c = (M/0)/(M/0) cr . Collapsing cores must 
have J c > 1 and ji c > 1, while cores with J c > 1 and 
/i c < 1 can become gravitationally bound but remain in a 
stable magnetostatic state (under ideal MHD conditions). 
Finally, cores with J c < 1 are Jeans-stable, and must re- 
expand after a transient turbulent compression, regardless 
of their value of /i c . Note that both J c and fi c are in gen- 
eral time-dependent, and also depend on the definition of 
the core's boundary. A fundamental question is then what 
kind of cores are produced in turbulent clouds as a func- 
tion of the turbulent parameters and of the cloud's own 
sub- or supercritical nature. 

In order to obtain insight on this, it is instructive to 
consider the formation of cores starting from idealized, ini- 
tially uniform-density clouds, truncated at a certain size, 
so that the cloud's mass-to-flux ratio is well defined. This 
represents also the periodic-boundary boxes of the numer- 
ical simulations. In reality, clouds are embedded in the 
diffuse atomic medium, and in principle there is always a 
large enough mass reservoir for them to be supercritical 
(e.g., Hartmann et al. 2001), but if the region size that 
needs to be considered is very large, then the time scale 
for accumulating this material in the molecular region is 
long, during which the cloud can be considered subcriti- 
cal over the evolutionary time scales of individual clumps 
within it. We thus consider both sub and supercritical 
uniform initial conditions, under ideal MHD, and a fixed 
global mass-to-flux ratio. 

In ideal MHD, it is clear that an initially subcritical core 
can become supercritical only if its parent cloud is super- 
critical. Indeed, the mass-to-flux ratio of an isolated cloud 
or core is independent of its radius, because both the mass 
and the flux are conserved quantities in this case. So, even 
if the whole cloud were to be compressed by a turbulent 
fluctuation, it could not increase its mass-to-flux ratio. As 
a consequence, a hard upper limit for the mass-to-flux ra- 
tio of a core is that of its parent cloud. Hence the need for 
the whole cloud to be supercritical if supercritical cores 
can ever arise under ideal MHD. 

The above hard limit on the mass-to-flux ratio of a core 
can also be seen as follows. Consider a cubic cloud of size 
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L with uniform density no and magnetic field Bq, and any 
flux tube of square cross section I 2 , with I <C L. Since in 
this case M/<f> = N /B , where No = n a L is the column 
density across the cubic domain, it is clear that M/<j> for a 
flux tube is independent of its cross section, and equal to 
that of the whole (uniform) cloud. If all the mass in the 
flux tube were to end up in a single core, the core would 
still be limited to having the same value of M/(f> as its par- 
ent, uniform cloud. If instead some fragmentation process 
(e.g., turbulent or gravitational) causes the formation of 
several cores along the flux tube, then each one of them 
will have a smaller value of M/(p than the parent cloud, 
because the sum of their masses is bounded from above by 
the mass of the flux tube, while the flux along the tube is 
constant. 

On the other hand, a lower estimate of the mass-to-flux 
ratio of a core can be obtained by considering simply a 
region of volume ^ 3 within the same uniform cloud. In this 
case, the mass-to-flux ratio of this parcel (which is not 
properly a core, since it is at the same uniform density 
as the whole cloud) scales as l/L. So, for a core of size I 
along the flux tube, we may expect its mass-to-flux ratio 
jLt c to lie in the range (1/L)/j,q < fi c < when its density 
n c lies in the range no < n c < (L/^) 3 no, where fio and no 
are respectively the parent cloud's mass-to-flux ratio and 
mean density. 

4. NUMERICAL EXPERIMENTS 

4.1. Method and parameters 

We use a total variation diminishing (TVD) scheme for 
solving the isothermal MHD equations in three dimen- 
sions. TVD is a second-order-accurate upwind scheme, 
and its implementation for isothermal flows is described 
in detail by Kim et al. (1999). Here we use an extension of 
this code including self-gravity and a random turbulence 
driver. The boundary conditions in all our simulations are 
periodic. We solve the Poisson equation using the usual 
Fourier method. In order to achieve second-order accuracy 
in time, we implement an update step of the momentum 
density due to the gravitational force as in Truclove et 
al. (1997). For the turbulence random driver, we follow 
the method in Stone, Ostriker & Gammie (1998). As dis- 
cussed in §2.2, we drive the turbulence at large scales, at 
wavenumber k — 2(2n/L), where L is the one-dimensional 
size of the computational box. We adjust the kinetic en- 
ergy input rate in order to maintain a roughly constant 
specified rms sonic Mach number M s . 

The simulations presented in this section can be con- 
sidered extensions of previous works by Pouquet, Leorat 
& Passot (1991), Ostriker et al. (1999), and Heitsch et al. 
(2001). In those studies, numerical simulations in which 
the entire computational domain is supercritical systemat- 
ically undergo collapse, albeit somewhat more slowly than 
non-magnetic cases (by factors of 2-3). Moreover, Li et 
al. (2004) have shown that all cores arising in their su- 
percritical simulations are supercritical by at least an or- 
der of magnitude. Only when the entire computational 
box is constrained to be subcritical is the collapse of both 
the large and the small scales prevented, giving rise to 
flattened structures. We extend those numerical works 
to study the lifetimes and criticality of the individual 
cores in relation to the criticality of their parent clouds. 



To this end, we now consider a series of simulations of 
Jeans-unstable numerical boxes, considering both magnet- 
ically sub- and supercritical cases. In view of our neglect 
of AD, AD- mediated contraction and subsequent gravi- 
tational collapse cannot occur in our simulations. How- 
ever, we can investigate whether magnetostatic, subcriti- 
cal cores with J c > 1 and yu c < 1 (cf. §3.3) form in the 
simulations, that could then evolve quasi-statically under 
AD. This is done in §5.2.1. 

The simulations are scale-free, and are characterized by 
three non-dimensional numbers: M s = a/c (the rms sonic 
Mach number, where a is the turbulent velocity dispersion 
and c is the sound speed), J = L/Lj (the Jeans number, 
giving the size of the box in units of the Jeans length 
Lj), and j3 = Pth/-Pmag = &^Poc 2 /B 2 (the plasma beta, 
defined as the ratio of the thermal to the magnetic pres- 
sures, where the subindex "0" denotes mean values in the 
box and c is the sound speed). The values of these param- 
eters for all runs are respectively given in columns 2, 3 and 
4 of Table 1, together with the run's resolution (column 
5) and suitable physical units, obtained as follows: The 
sound speed is assumed to be fixed, at c = 0.2 km s _1 
(implying T = 11.4 km s _1 ). Then, a Larson- like relation 

of the form M s = 5 (L/lpc) 1/2 (e.g., Larson 1981; Blitz 
1993) gives the box size (column 6) as a function of M s as 

L = 0.04M 2 pc. (5) 

The Jeans number and the box size give the mean number 
density density (column 7) as 

2 



n = 500 



cm 



(6) 



_(L/l pc)_ 

where we have assumed a mean molecular weight m = 
2Amn- In turn, the mean density and j3 (at the above 
value of c) give the mean field strength (column 8) as 

Bo = 0.205 ' M G. (7) 

The last two columns in Table 1 give the rms turbulent 
velocity dispersion and the mass-to-flux ratio of the simu- 
lation (in units of the critical value), /i = {13/ f3 c ^) 1 ^ 2 , where 
P cr is given as a function of J as 

/3 cr ~ 2/ (tt 2 J 2 ). (8) 

This in turn is obtained assuming that the critical mass-to- 
flux ratio is (M/<p) CI w (4tt 2 G)~ 1/2 (Nakano & Nakamura 
1978). Specifically, we get /3 cr w 0.023 for J = 3 and 
/3 cr w 0.013 for J = 4. We have verified numerically the 
critical (3 value for the J = 3 case as follows. We set up a 
white-noise initial velocity field with rms amplitude 10 _4 c 
in each velocity component, and observe whether the max- 
imum density in the simulation increases or becomes satu- 
rated in time. Indeed, the cases with (3 > 0.03 collapse in 
a runaway fashion, whereas the cases with (3 < 0.01 form 
magnetically-supported configurations, flattened along the 
mean field direction. 

4.2. Resolution considerations 

It is important to know when our results can be con- 
sidered robust physical ones rather than numerical arti- 
facts caused by insufficient resolution. A useful criterion 
for guaranteeing that the gravitational collapse of an ob- 
ject is followed adequately was given by Truelove et al. 
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(1997). This criterion requires that the minimum Jeans 
length within the simulation be resolved by at least four 
zones. Heitsch et al. (2001) increased this requirement 
to six zones for the proper resolution of MHD waves, al- 
though this may be an excessive constraint, as those au- 
thors showed that the magnetic field can only prevent 
collapse if it is strong enough to provide magnetostatic 
support. Our interest just being to determine whether 
magnetostatically-supported structures arise, we consider 
that fulfillment of the 4-zone Jeans condition is sufficient. 
This is consistent with the results described below. 

In general, the constraint of resolving the Jeans length 
by Nl 3 zones limits the maximum resolvable density to 

2 

(9) 



77, 



JN Lj 



no, 



where A^ox is the number of zones per dimension in the 
computational box. Note that our definition of J is dif- 
ferent from that of Truelove et al. (1997). For the Jeans 
criterion, Nl., — 4, so for our J = 4 runs, this gives n rcs = 
256no. For the more stringent criterion of MHD-wave res- 
olution, Nlj = 6 zones, in which 114n . 

4.3. Core analysis procedure 

In the §5.2, we describe in detail the evolution of the 
magnetic criticality and Jeans number of some selected 
cores in the simulations. Here we describe the proce- 
dure for measuring their mass-to-flux ratio p c and Jeans 
number J c . Note that this is not a fully unambiguous 
task, because the medium is a continuum. Thus, the 
cores' boundaries, and therefore their masses and sizes, 
are somewhat arbitrary. Moreover, the cores are in gen- 
eral part of a larger parent clump, which is also involved 
in the dynamics. A precise measurement of the struc- 
tures' energy balance would require to measure in prin- 
ciple all the terms in the virial theorem for them, in- 
cluding the surface terms (Ballesteros-Paredes & Vazquez- 
Semadeni 1997; Ballesteros-Paredes et al. 1999a; Shad- 
mchri, Vazquez-Semadeni & Ballesteros-Paredes 2001; 
Tillcy & Pudritz 2004). Thus, here we only measure the 
mass-to-flux ratio of the cores in an approximate way, as 
follows. 

First, we define a structure (clump or core, which here 
we will generically refer to as "core") as a connected re- 
gion around a local density peak, in which the density is 
above a threshold value n t . Then, the mass of the core is 
computed directly from the density field as M c = J v pdV, 
where V c is the core volume. To compute the magnetic 
flux within the core, we need in principle to first compute 
the (vector) average of the magnetic field ((B)) within the 
core's volume, and then define a surface S normal to (B) 
that bisects the core, in order to compute C = f s B • hdA, 
where n is the unit vector normal to S. For our purposes 
here, we use a simpler, approximate procedure. For su- 
percritical runs, in which the cores are not very far from 
round, we define the core radius as R c = (Wc/Air) 1 / 3 . It 
is easy to show that, if instead of spherical, the core is a 
cylinder of volume irR 2 dz, then an aspect ratio a = R/dz 
gives an error of a factor of a 1 / 3 in the above estimate 
of the radius. At most, the aspect ratios of the cores in 
our moderately supercritical run are < 4, and are even 
smaller in the strongly supercritical run. Thus, the errors 



in the estimate of R are by factors < 1.6. For the sub- 
critical run, and very flattened cores, we define the core 
radius as the maximum distance between the density peak 
and all points belonging to the core. Typical sizes for the 
"clumps" (defined by n t = 10no) and for the "cores" (de- 
fined by n t — 10n ) are ~ 0.08 and 0.04 times the box size 
(~ 0.32 and 0.16 pc). 

With these definitions, we finally estimate the flux 
through the core is <p c w ttBR 2 , where B is the mean 
magnetic field intensity through the core. The core's 
Jeans number is simply computed as J c = i?/Lj. c , where 

1/2 

Lj, c = (ttc 2 /G/9 c ) , is the mean Jeans length in the core, 
and p c is the mean density in the core. 

5. NUMERICAL RESULTS 
5.1. Global analysis 

Figure 2 shows the evolution of the global density maxi- 
mum n max for the four simulations in Table 1 with M s = 10 
and J = 4 and (3 = 0.01, 0.1, 1 and oo, respectively labeled 
M10J4/3.01, M10J4/3.1, M10J4/31 and M10J4/?oo. The 
time axis is shown both in Myr (lower axis) and in units 
of the global free-fall time (if g = Lj/c, upper axis). Note 
that runs M10J4/3.01 and M10J4/3.1 are evolved for over 
10 Myr. This is a long time compared to recent estimates 
of cloud lifetimes of a few Myr (Ballesteros-Paredes, Hart- 
mann & Vazquez-Semadeni 1999; Hartmann et al. 2001). 
The long integration times are just for statistical purposes, 
and by no means imply that we advocate longer lifetimes. 

From fig. 2 it is seen that in the subcritical ([3 = 
0.01) case, n max is typically ~ 40no and seldom exceeds 
~ 100no. The same behavior is observed in the other 
subcritical runs (not shown), M10J3/3.01, M8J4/3.01, and 
M10J4/3.01LR. Thus, in these simulations core-like densi- 
ties (n > 200no ~ 10 5 cm~ 3 ) arc not produced. The fact 
that the maximum density does not stay constant also sug- 
gests that there are no magnetostatic clumps. Note that 
these results cannot be attributed to numerical diffusion, 
whose effect would be expected to act in the opposite di- 
rection of causing artificial flux loss and possibly collapse. 

In contrast, fig. 2 also shows that in the supercritical 
runs n max reaches values ~ 10 4 no after roughly 1/2 tf e , 
although the rise from densities > 100n to ~ 10 4 n oc- 
curs on a much shorter time scale. This corresponds to the 
first local collapse event, as discussed in §4.3. Finally, the 
non-magnetic run M10J4/3oo is seen to start producing col- 
lapsed objects much earlier than the magnetic supercrit- 
ical runs. Thus, both the supercritical runs and the non- 
magnetic ones produce the first collapsed objects on times 
shorter than the global free-fall time as a consequence of 
the production of locally gravitationally unstable objects 
by the turbulence. However, it is seen that the effect of 
the turbulent compressions is diminished by the presence 
of the magnetic field even in the strongly supercritical case 
in comparison to the non-magnetic one. 

5.2. Evolution of core properties 

5.2.1. Subcritical case 

In fig. 3, we show an animation of the evolution of 
the subcritical run M10J4/3.01. Depicted are three iso- 
surfaces at densities 10 (blue), 30 (yellow) and 100 (red) 
times n . The time interval between successive frames is 



Lifetimes and Evolution of Molecular Cloud Cores 



9 



At = 0.002 t s , where t s = 20 Myr is the sound crossing 
time for the M s = 10 (L = 4 pc) boxes. In this simulation, 
most density peaks ( "clumps" , in this case) are seen to be 
extremely transient in nature, with lifetimes 0.01-0.075 t s , 
or 0.2-1.5 Myr (some 5-30 frames), and with poorly de- 
fined identities, the clumps actually merging with other 
structures or simply morphing into different shapes and 
sizes. 

In this simulation, we only observe one clump that be- 
comes gravitationally bound and survives for a few local 
free-fall times. It appears near the left boundary of the 
box around frame 205, and moves across it to remain near 
the right boundary for the rest of its existence. This core 
maintains peak densities n pe ak > 100no = 5 x 10 4 cm~ 3 
during 7 frames, from frame 221 to 227 (0.28 Myr), and 
ftpeak ^ 60no during <~ 28 frames (~ 1.1 Myr), from frame 
210 to frame 237. Figure 4 shows the evolution of /z c and 
J c for this clump (defined by n t = 40no) from frame 210 to 
240 (t = 8.4 to 9.6 Myr). We see that it remains strongly 
subcritical (well below the value for the box of 0.9), but 
it manages to remain significantly super-Jeans for slightly 
over 1.2 Myr. 

This clump then represents the closest our subcritical 
simulation gets to producing the conditions for subsequent 
quasi-magnetostatic evolution mediated by AD. Neverthe- 
less, we see that this clump's lifetime is still quite short, be- 
cause it is not very strongly gravitationally bound, and the 
turbulence manages to tear it appart after having formed 
it. It is then natural to ask whether AD would have enough 
time, had it been included in the simulations, to allow it 
to become more strongly bound, and ultimately collapse, 
before it is dispersed again by the turbulent motions. 

To answer this question, Appendix I presents a simple 
estimate of the classical AD time scale under the condi- 
tions of this clump. It is found that tad ^1-3 Myr, which 
is very close to the duration of this clump in our simula- 
tion, suggesting that AD is possibly capable of "capturing" 
it and helping it decouple from the external turbulence. 
Note, however, that this value of tad is not significantly 
larger than the collapse times we observe in the supercrit- 
ical simulations, and thus a core that evolves on this time 
scale can hardly be considered quasi-static. 

Another relevant question is whether such clumps would 
approach a magnetostatic state if the turbulence were de- 
caying. To test this possibility, we have conducted the 
additional experiment of turning off the turbulent driving 
at frame 210, right after the clump forms. In this case, we 
find that indeed the clump undergoes a few bounces, but 
maintains its identity for at least 1.6 Myr (the duration 
of the integration we performed), with no sign of disrup- 
tion towards the end. This suggest that, if turbulence is 
turned off right after it forms the clump, then the clump 
does evolve towards a magnetostatic state, and can easily 
be captured by AD, to evolve in accordance to the stan- 
dard model. Presumably, other clumps might undergo the 
same evolution in a decaying situation, as in the simula- 
tions by Li & Nakamura 2004. 

We conclude from this section that even the densest, 
longest-lived clumps arising in our subcritical simulation 
are dispersed here because the strong magnetic field does 
not allow them to collapse nor achieve very strong grav- 
itational binding, and makes them susceptible of subse- 



quent disruption by the same turbulent field. The longest 
clump lifetimes are marginally long enough for them to 
undergo AD-mcdiated evolution, albeit on relatively short 
time scales. However, such clumps could easily become 
magnetostatic in a decaying turbulent field. The funda- 
mental question then becomes whether molecular clouds 
are continually driven or decaying (see discussion in §6). 

5.2.2. Moderately supercritical case 

Figure 5a shows the evolution of the moderately super- 
critical run M10J4/3.1, again with At = 0.002 t s , but now 
with the isodensity contours denoting 10 (blue), 100 (yel- 
low) and 1000 (red) xn . It is seen that in this run, objects 
reaching much higher peak densities, n poa k > 5000no, form 
since early in the evolution. Specifically, three such objects 
form at times ~ 1.76, 5.2 and 8.4 Myr (frames 44, 130 and 
210, respectively). Once they have reached those densi- 
ties, these objects are already collapsed, as they simulta- 
neously have J c ^> 1 and /i c 3> 1, and so nothing sup- 
ports them against their self-gravity. In particular, they 
have (U c > ^box, indicating that numerical diffussion has 
caused magnetic flux loss within them. They do not pro- 
ceed to a singularity because, once most of the collapsing 
mass is in a single, completely unresolved numerical cell, 
the collapse cannot be followed any further by the numer- 
ical grid. These cells, which we refer to as "sink cells", 
can be thought of as the fixed-grid equivalent of sink par- 
ticles in SPH simulations (Bate, Bonnell & Price 1995), 
and are standard fare in fixed-grid numerical studies (e.g., 
Hcitsch ct al. 2001; Li ct al. 2004). Note that, however, 
the sink cells continue to accrete from their surroundings 
at a slower rate, and eventually the simulation finally stops 
when the density gradient becomes too large. Neverthe- 
less, as we show below, the onset of collapse occurs under 
well-resolved conditions. 

Now we focus on the evolution of the first collapsing 
object. From the animation in fig. 5a, we see that the 
collapsing core actually forms from the merger of two pre- 
vious cores, which had formed almost simultaneously at 
frame 36 within a clump located near the center of the 
computational box. The two cores merge at frame 44, 
which marks the "birth" of the collapsing core, all hap- 
pening within the larger clump. Thus, we investigate the 
criticality and stability of both the "parent" clump and the 
"daughter" cores at frames 30, 40 and 50. Figures 6a and 
6b respectively show 3D iso-density surface maps of the 
clump/core system in question at frames 30 (t = 1.2 Myr) 
and 40 (t = 1.6 Myr). The iso-surfaces are drawn at den- 
sities n = 10n (blue), n — 40n (yellow) and n = 100n 
(red). 

At frame 30 (fig. 6a), it is seen that the n = 100no iso- 
surface does not exist, since the peak density is < 100n . 
It is also seen that the n = lOrio isosurface defines a sin- 
gle structure, which we refer to as the parent "clump", 
while the n = 40no isosurface resolves the two cores. At 
frame 40 (fig. 6b), the n = 40n isosurface does not resolve 
the two daughter cores anymore, although the n = 100n 
isosurface still does. It is important to remark that the iso- 
surfaces drawn at a certain value of nt at different times 
do not trace the same gas parcels; i.e., the isodensity sur- 
faces are not Lagrangian boundaries. Neither are clump 
boundaries defined by means of standard clump-finding 
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algorithms. 

Figure 7 shows the evolution of the mass-to-flux ratio 
(panel a) and of the Jeans number (panel b) for this col- 
lapsing object from frame 30 to frame 50 (t = 1.2 to 2 
Myr). Three curves are shown in each panel, each corre- 
sponding to different values of the clump-defining thresh- 
old density n t . Note that two cores exist at frame 30 for 
n t = 40n , but their fi c and J c are nearly identical, and 
so they are described by the same (dotted) line. Similarly, 
at frame 40 there are two cores defined by the nt — lOOrto 
level, and in this case their /i c and J c are different, so 
they are described by two separate (dashed) lines, which 
in fact start at frame 40, because the cores did not exist 
at threshold nt = lOOrio at frame 30. 

From fig. 7 it is seen that that the parent clump is already 
Jeans-unstable and magnetically supercritical at frame 30, 
even though it contains two well-defined substructures, 
both of which individually arc subcritical and sub- Jeans at 
this time, and are therefore probably not collapsing them- 
selves at this time. In this sense, the merger of the two 
cores can be considered a part of the global collapse of 
the larger structure, a fact which is confirmed by noting 
that in the animation, the size of the n — lOrto isosurface 
is also shrinking, implying that most of the mass in the 
whole clump is involved in the collapse (see also fig. 6). 

Note that the super- Jeans and supercritical character of 
the n t = 10no clump at frame 30 cannot be attributed to 
numerical diffusion, as its mean density is only ~ 20rio, 
with its peak only reaching 83no, thus very safely below 
even the most stringent (and probably excessive) resolu- 
tion constraint of Heitsch et al. (2001) of n rcs = 114no. 
This is also supported by the fact that at this time the 
mass-to-flux ratio of the clump, albeit supercritical, is still 
smaller than that of the whole box, in agreement with the 
constraint that no structure within the simulation can have 
a value of n larger than that of the whole box under ideal 
MHD. The latter property holds also for the two cores de- 
fined by n t = lOOrio at frame 40, whose mean and peak 
densities are respectively (ni) = 153rio, n pcli ^ 1 = 337no 
and (n 2 ) = 121n , n poa k 2 = 166n . It is only at frame 50, 
when the core has a peak density of over 5000no and an 
average density of 329n (above n t = 40n ), that its mass- 
to-flux ratio exceeds that of the box, indicating that at 
this time numerical diffusion has clearly caused significant 
loss of magnetic flux. However, the collapse was initiated 
when the whole structure was cleanly resolved, and there- 
fore must be considered a robust physical result. That 
is, the formation of a supercritical, super-Jeans structure 
causes collapse and ultimately numerical diffusion, rather 
than the collapse being an artifact of numerical diffusion. 
The same happens during the formation and collapse of the 
other two collapsed cores (at frames 130 and 210). This 
conclusion is also supported by the fact that no collapse 
occurs in the subcritical simulation, even though it is only 
marginally so (fi = 0.9), which shows that the simulations 
accurately do not allow collapse when it physically should 
not occur. 

Finally, two more points are worthy of notice. One is 
that this first collapsing object goes from first appear- 
ance to a fully collapsed state in less than 2tf c , where 
t{ c = Lj tC /c is the local free-fall time in the core, with 
Lj iC being the Jeans length at the core's mean density. At 



a mean density (n) c = lOOrio, tf c = 0.5Myr = 12.5 frames. 
The collapse corresponds to the rapid rise of n max from 
~ 10 2 to <~ 10 4 no in fig. 2. The second is that, as in the 
non-magnetic case, moderate-density (n poa k > 10 2 cm~ 3 ) 
chimps arise that redisperse on time scales 1/2-1 tf c (< 10 
frames). 

5.2.3. Strongly supercritical case 

In the strongly supercritical run M10J4/31 (fig. 5b), 
collapsing objects undergo a very similar evolution as 
in the moderately supercritical one, M10J4/3.1. How- 
ever, this run exhibits an additional interesting phe- 
nomenon: Besides the rapidly-collapsing dense cores and 
the rapidly-dispersing intermediate density clumps dis- 
cussed above, we also observe the formation of two 
long-lasting, intermediate-density clumps (n poa k 200-300 
cm~ 3 ) which do not collapse, but which take times ~ 
0.15 t s (~ 0.6 if g , ~ 75 frames, or ~ 3 Myr) to redisperse. 
These are seen as the long-lasting clumps that do not have 
a red isodensity surface in the animation. One of them is 
formed at frame 16 (t = 0.64 Myr) near the top of the 
simulation, and disperses around frame 84 (t = 3.36 Myr). 
The other one forms at frame 42 (t = 1.68 Myr) towards 
the left of the simulation, at roughly half the box height, 
and disperses around frame 120 (t = 4.80 Myr). These 
clumps are marginally resolved according to the Jeans cri- 
terion, but nevertheless end up redispersing, rather than 
collapsing, suggesting that numerical diffusion of the mag- 
netic field is not a concern yet in these clumps. 

These clumps last a time about twice the longest dura- 
tions seen in the subcritical case, and moreover have very 
well-defined identities, contrary to the situation in the sub- 
critical case. It is thus interesting to investigate their phys- 
ical properties. Figure 8 shows the evoulution of /j, c and 
J c for the first of these clumps from frame 20 to frame 70 
(t = 0.5 to t — 3. Myr). It can be seen that this core is su- 
percritical (/i > 1), but sub- Jeans (J c < 1) at all times and 
values of the core-defining threshold density n t , except at 
frame 20 and n t — 10n , at which J c is marginally larger 
than 1. Given the approximate nature of our calculations 
(recall the core's radius is defined assuming a spherical 
geometry), and that turbulent support is not included in 
our analysis, we consider that the value J c = 1.06 found at 
that particular time and density threshold is still consis- 
tent with the core not undergoing collapse in practice. We 
conclude that these somewhat longer-lived, non-collapsing 
cores come very close to collapse, but miss it, explaining 
their longer lifetimes in terms of our estimate of §3.2. It 
appears unlikely that the inclusion of AD in the calcula- 
tions could help these cores to collapse, as they are already 
supercritical, and so their support is mainly thermal. 

It is also interesting that these longest-lived cores arise 
in the most strongly supercritical simulation, in which 
the magnetic support is weakest. A possible explana- 
tion is that, somewhat countcrintuitivcly, turbulent flows 
with weak magnetic fields tend to develop more moderate 
density fluctuations than cither non-magnetic or strongly 
magnetized ones, as has been observed by various groups 
(Passot et al. 1995; Ostriker et al. 1999; Heitsch et al. 2001; 
Ballesteros-Paredes & Mac Low 2002). Thus, the produc- 
tion of moderate-amplitude density fluctuations that can 
"land" close to equilibrium without directly overshooting 
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to collapse may be more likely in the weak-field case. 

5.2.4. Non-magnetic simulation 

Figure 9 shows an animation of the nonmagnetic run 
M10J4/3oo. This simulation stops after only 45 frames, 
indicating that the collapse is more violent than in the 
magnetic cases. We can only give lower bounds to the 
clump lifetimes in this case, and a more detailed descrip- 
tion must await a more robust handling scheme for the 
collapsed objects, to be presented elsewhere. 

Two features are worth noting here. First, a significantly 
larger number of collapsed objects and clumps have formed 
by the end of the run than at the same time in either MHD 
run. This result must be considered preliminary only, as 
it is necessary to verify that the larger number of cores 
is a real effect, and not just a statistical fluctuation. For 
example, Heitsch et al. (2001) reported that the statistical 
fluctuations between realizations with identical parame- 
ters but different random initial conditions and driving 
were larger than the systematic differences between mag- 
netic and non-magnetic cases. However, they did not dis- 
cuss the number of cores formed in their simulations. The 
difference in core number between the magnetic and non- 
magnetic runs in our simulations is, however, striking, and 
we plan to investigate this in a forthcoming paper consider- 
ing mean trends with the global parameters for ensembles 
of simulations. 

The second point to note is that the time scales from 
formation to full collapse seem to cover a wider range, of 
<~ 1/2 to 1 Myr, with the fastest-collapsing objects being 
those that involve collisions of clumps. Moreover, a pop- 
ulation of non-collapsed objects is seen by the end of the 
run for which we cannot observe the final fate nor the cor- 
responding time scale. Several re-expanding clumps are 
also observed. The investigation of these trends must also 
await a more robust handling scheme for the sink cells. 

Nevertheless, these preliminary results lead us to spec- 
ulate that the presence of a magnetic field can further re- 
duce the efficiency of collapsed-object formation in com- 
parison with that of an equivalent non-magnetic case 
mostly by reducing the probability of core formation, 
rather than by significantly delaying the collapse of in- 
dividual objects. 

6. DISCUSSION, IMPLICATIONS AND CAVEATS 

Several interesting points can be noted about the results 
described above: 

1. Our results suggest that, at least under the assump- 
tions of our numerical simulations, molecular cloud clumps 
and cores are in general likely to be dynamic, out-of- 
cquilibrium structures, rather than quasi- hydro/magneto- 
static structures. This is however not in conflict with 
observations, as several studies show that synthetic ob- 
servations of the clouds and the clumps within them in 
numerical simulations compare favorably to observations 
(e.g., Padoan et al. 1999, 2003; Klessen 2001; Ostriker et 
al. 2001; Ballesteros-Paredes & Mac Low 2002; Gammie et 
al. 2003; Klessen, Ballesteros-Paredes, Vazqucz-Scmadcni 
& Duran 2004; Schmeja & Klessen 2004). In particu- 
lar, dynamically- formed cores can have angle- averaged col- 
umn density profiles that resemble that of a BE-sphere 
(Ballesteros-Paredes, Klessen & Vazquez-Semadeni 2003), 



or in fact various other functional forms (Harvey et al. 
2003). This is because much detail is lost by the line-of- 
sight and angle averaging. 

2. The clump/core lifetimes that we have found in this 
paper (1/2-1 Myr for collapsing objects, 1/2-3 Myr for 
non-collapsing ones) also compare favorably with obser- 
vational estimates of core lifetimes, such as those of Lee 
& Myers (1999), who obtained estimated lifetimes of 0.3- 
1.6 Myr. However, an interesting remark is in order. Lee 
& Myers obtained those estimates under the assumption 
that all cores undergo AD-mediated contraction, infall, 
and then protostellar emission, so that the number of star- 
less cores (tisl) and of cores with young stellar objects 
(YSOs) (riYSo) are respectively proportional to the time 
spent in the prestellar and protostellar epochs during this 
process. However, our results advocate the possibility that 
not all clumps undergo collapse, and in this case it is easy 
to show that the ratio hsl/^yso is an upper bound to the 
lifetime ratio. Indeed, define e = tisl/^yso- If the num- 
ber of starless cores comprises a fraction e of failed cores 
and a fraction 1 — e of true pre-stellar cores (with number 
rips), then 



nyso nyso \1 



(10) 



Now, since the ratio of true pre-stellar objects to proto- 
stellar ones is indeed equal to the ratio of pre-stellar to 
protostellar lifetimes tps/tyso, we find 



TPS 



tyso \1 - e 



1 



> 



TPS 
TYSO 



(11) 



which is the desired result. Although the actual fraction 
of failed cores has not been quantified in our numerical 
simulations (again, ensembles of simulations are needed) , 
and is unknown in the observations, this consideration is 
consistent with the fact that the longest lifetimes of col- 
lapsing objects in our simulations are shorter by a factor 
of - 2/3 than the upper bound of Lee & Myers' (1999) 
estimates. 

3. It is important to compare the fraction of collaps- 
ing versus non-colapsing clumps in the simulations with 
observational data. This will be attempted in detail in a 
separate paper, but here we can give the following sim- 
ple estimate. Assuming that a minimum number density 
of <~ 10 4 cm~ 3 is necessary to excite a high-density tracer 
such as ammonia, a typical "ammonia core" consists of the 
"tip of the iceberg" region of a clump above this minimum 
density. Such minimum density can be identified with 
our core-defining threshold density n t , giving n t ~ 20n . 
In our simulations, typical failed, re-expanding cores de- 
fined with this value of nt only reach mean densities of 
up to 2-3 times n t (~ 50n = 2.5 x 10 4 cm" 3 ). Thus, 
it seems that most of the re-expanding clumps and cores 
in our simulations would lie below the ammonia detection 
limit, and would preferentially appear in lower-density- 
tracer surveys, while surveys in high-density tracers such 
as ammonia may tend to preferentially select collapsing 
objects. Detailed synthetic observations of the simulations 
with various clump-defining density thresholds are neces- 
sary to determine whether the observed fraction of inward 
versus outwards (i.e., contracting vs. re-expanding) mo- 
tions in cores (e.g., Lee, Myers & Tafalla 2001) are well 
reproduced by the simulations. 
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4. In the analysis of §5.2, we have neglected turbu- 
lent support of the clumps and cores. Interestingly, the 
{J c ,Hc) parameter space seems to describe very well the 
outcome of the clump/core system evolution. This may 
be because at the small sizes of these objects, the internal 
turbulence may be of secondary importance. Indeed, real 
cores are typically transonic or even subsonic, a feature 
reproduced by the simulations (Klessen et al. 2004). In 
this case, the turbulent support is comparable or smaller 
than the thermal support, and so only the magnetic and 
thermal supports remain for consideration. 

5. The caveat that our results hold only within the realm 
of the assumptions described in §2 should be stressed. Al- 
though in that section we have referred to observational 
evidence suggestive that the numerical setup and parame- 
ters used in this paper are the most realistic (except for the 
neglect of AD), the issue of whether molecular clouds are 
nearly isothermal, continually driven and supercritical is 
still unsettled. Our clump lifetimes in the subcritical case 
are short because the turbulence is continually driven. We 
cannot overemphasize the importance of continued obser- 
vational surveys directed at determining the true condi- 
tions in molecular clouds, or the fraction of clouds under 
each type of condition. Also, numerical simulations of the 
interstellar medium at large, in which the formation of 
molecular clouds is accurately simulated, can shed light 
on this issue. 

7. SUMMARY AND CONCLUSIONS 

In this paper we have discussed the lifetimes of clumps 
and cores in both the non-magnetic and magnetic cases, 
in single-phase, isothermal, turbulent models of molecular 
clouds. Under these conditions, we have remarked that, 
in the non-magnetic case, the hotter, low-density gas nec- 
essary for the confinement of stable BE-type hydrostatic 
configurations is absent, and so we do not expect these 
structures to form out of transient turbulent fluctuations 
in the clouds. We have then given a qualitative discus- 
sion suggesting that non-truncated ("extended") isother- 
mal structures may be expected to be unstable in general, 
in agreement with the fact that non-magnetic polytropic 
simulations with 7 e < 4/3 do not form hydrostatic struc- 
tures (e.g., Vazquez-Semadeni et al. 1996; Klessen et al. 
2000; Bate et al. 2003; Li et al. 2003). Thus, we have 
suggested that in the non-magnetic case, clumps should 
be expected to be completely transient structures, either 
proceeding to collapse right away upon the turbulent com- 
pression that formed them, or else re-expanding to merge 
back into their environment ("failed" clumps). 

For failed clumps, wc have given a simple estimate for 
the characteristic re-expansion time, as slowed down by 
the self-gravity of the clump, finding times ~ a few clump 
free-fall times, comparable to the lifetime estimates of Lee 
& Myers (1999). The longest lifetimes correspond to re- 
expanding clumps that come very close to equilibrium be- 
tween self-gravity and thermal pressure. In general, failed 
clumps are not destined to form stars. However, in the 
simulations there is always a population of these clumps 
present, in agreement with some observational suggestions 
(Taylor et al. 1996). These can be associated with at least 
part of the observed starless cores, which traditionally have 
been thought of as being in a pre-stellar phase. 



In the magnetic case, we have presented a set of of 
drivcn-turbulcncc, idcal-MHD numerical simulations fo- 
cusing on the lifetimes and criticality of the individual 
clumps and cores in relation to the conditions in their par- 
ent clouds. We have found that: 

a) Magnetostatic clumps do not form in the case of sub- 
critical environments; the clumps that form in this case 
have mean densities n <~ 10-100n (5 x 10 3 -5 x 10 4 cm -3 ), 
and lifetimes ~ 0.2-1.3 Myr, ultimately being dispersed 
by the turbulence in the environment. A few of these 
clumps nevertheless become moderately gravitationally 
bound, being the ones with the longest lifetimes. An esti- 
mate of the AD time scale in these clumps gives character- 
istic times > 1.3 Myr, suggesting that they can marginally 
be "captured" by AD and increase their mass-to-flux ratio 
before they can be dispersed, thus increasing their gravita- 
tional binding, to eventually become supercritical and col- 
lapse. Thus, they adjust to the standard model scenario, 
except for the fact that with these time scales, AD does 
not significantly increase the cores' lifetimes in comparison 
to those of supercritical or non-magnetic ones. This ques- 
tions the usefulness of AD-mediated contraction as a SFE- 
decreasing agent in the standard model, and the criticality 
of the parent may become irrelevant for the formation of 
supercritical cores in short time scales. 

b) In the supercritical case, dense cores (n ~ 100-1000 
xrto, or 5 x 10 5 -5 x 10 6 cm~ 3 do form, being supercrit- 
ical since their early stages and collapsing on time scales 
<~ 2if c <~ 1 Myr. Ambipolar diffusion is not required for 
collapse in this the objects rapidly become super- 
critical by the turbulent compressions themselves, since 
there is a large enough mass reservoir around them. The 
lifetimes in the magnetic case are in general again in agree- 
ment with the observational estimates of Lee & Myers 
(1999) and our own simple, non-magnetic estimate. 

c) As in the non-magnetic case, coexisting with the 
collapsing cores, numerous lower-density (n < lOOrio) 
"failed" ones form and disperse in times comparable to 
their free-fall times. These "failed" clumps include the 
longest-lived objects observed in our simulations, again 
in agreement with our simple calculation, which suggests 
that the longest lifetimes correspond to re-expansion from 
states very close to equilibrium. However, it appears un- 
likely that the inclusion of ambipolar diffusion could in- 
duce their collapse, as these objects are already supercrit- 
ical, and the reason they do not collapse is because they 
do not manage to become gravitationally unstable with 
respect to the thermal support. 

Within the framework of our assumptions, our results 
support the developing notion that star formation is a dy- 
namic process, in which the only hydrostatic configura- 
tions may be the stars themselves. 
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Lifetimes and Evolution of Molecular Cloud Cores 



Table 1 

Numerical simulation parameters. a 



Run name 


M s 


J 


a 


Resolution 


L(pc) 


n (cm 3 ) 




er(km s x ) 


b ^ 


M10J4/3oo 


10 


4 


oo 


256 3 


4 


500 


0.0 


2.0 


oo 


M10J4/31 


10 


4 


1 


256 3 


4 


500 


4.6 


2.0 


8.8 


M10J4/3.1 


10 


4 


0.1 


256 3 


4 


500 


14.5 


2.0 


2.8 


M10J4/3.01 


10 


4 


0.01 


256 3 


4 


500 


45.8 


2.0 


0.9 


M10J4/3.01LR 


10 


4 


0.01 


128 3 


4 


500 


45.8 


2.0 


0.9 


M10J3/3.01 


10 


3 


0.01 


128 3 


4 


281 


34.4 


2.0 


0.7 


M8J4/3.1 


8 


4 


0.1 


128 3 


2.56 


1220 


22.6 


1.6 


2.8 


M8J4/3.01 


8 


4 


0.01 


128 3 


2.56 


1220 


71.6 


1.6 


0.9 



a See description in §4.1. 

b Turbulcnt velocity dispersion in the computational box. 
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APPENDIX 

AMBIPOLAR DIFFUSION TIME SCALE IN A SUBCRITICAL CLUMP 

Here we estimate the ambipolar diffusion (AD) time scaie tad under the conditions of the longest-lived clump in the 
subcritical run M10J4/3.01. To this end, we use the formulas given by Fiedler & Mouschovias (1992). They write 

8 4 f 8 riy 1 

where r s = R c /c s is the sound crossing time and rg is free-fall time scale 

37T 



T «=\m^) ■ (A2) 

Note that this definition of the free-fall time is smaller by a factor of y/3/32 than that used in the main body of this 
paper. This, however, causes no inconsistency, as we compare the AD time scale and core lifetimes in physical units. 
Finally, r nl is the neutral-ion collision time scale, given by 

r ni = (A3) 

in which m n and m\ are the neutral and ion particle masses, respectively, p n and p\ the neutral and ion mass densities, 
and (aw) i n is the neutral-ion collision rate. We can write the above equation as 

1+ rnn 

rm = -r^t- (A4) 
Fiedler & Mouschovias also give an approximate expression for the ion number density as 

. 10 5 cnr 

with canonical values K = 3 x 10~ 3 cm -3 , k = 1/2, and (aw)i n — 1.69 x 10 _9 cm 3 s _1 . We also take m\ ~ 29.0toh and 
r7i n = 2.36mn, with mn = 1.67 x 10~ 24 g. 

Defined out to a threshold density n t = 40no, the core in question has a mean density of 55rio = 2.75 x 10 4 cm -3 at the 
time when it attains the maximum vsIug of its m£tss-to-flux ratio, /Xc,max = 0.28 (see fig. 4). This core is highly flattened, 
and thus we compute its radius as the maximum distance between the density maximum and the clump's periphery (cf. 
§4.3), finding R c = 0.3 pc. With these values we thus obtain r ni = 1.29 x 10 4 yr, t b = 2.02 x 10 5 yr, r s = 1.5 x 10 6 yr, 

and [l - (8/7r 2 )(T ff /T s ) 2 ]~ = 1.015. Substituting in eq. (Al) we finally obtain 

tad ps 9.2-r ff = 1.9 Myr. (A6) 
This estimate is very close to the simpler, more general estimate often quoted (e.g., Ciolek & Basu 2001) of tad ~ 
T 2 /(1.4r ni ) = 11.2T ff = 2.3 Myr. 

A correction to this estimate to take into account the initial value of the mass-to-flux ratio in the core, /j c has recently 
been given by Ciolek & Basu (2001). From their fig. 1 and the value k = 1/2 that we have used here, the correction factor 
is C(/i c o) ~ 0.7. This is the smallest possible value of their correction factor at this value of fj, c . So we conclude that 

tad > 1-3 Myr. (A7) 
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1.0 1.2 1.4 1.6 1.8 2.0 2.2 



Fig. 1. — a) (Left panel) Re-expansion time of a core, in units of the free-fall time, denned as the time necessary to double its initial 
radius, as a function of the initial radius r%, normalized to the equilibrium radius, b) (Right panel) Simulated "observer's perspective" of the 
re-expansion time, giving the time to go from the initial radius r\ to a final one with density 0.1 p c , where p c is the density of the equilibrium 
configuration. This emulates the time the core would be visible to a high-density tracer sensitive to densities > 0.1p c . 
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Fig. 2. — Time-evolution of the global maximum of the density field for runs M10J4/3.01 (dotted line), M10J4/3.1 (solid line), M10J4/31 
(short-dashed line) and M10J4/3oo (long-dashed line). The numbers indicate the corresponding values of fi. The upper time-axis is given in 
units of the global free-fall time of the simulation. 
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Fig. 3. — Animation of the density field for the subcritical run M10J4/3.01, showing transparent density iso-surfaces at n = lOrao (blue), 
n = 30no (yellow), and n = lOOnrj (red). The figure shown in the paper edition corresponds to frame 223, at which the moderately-long-lived 
core attains its time- wise maximum central density, of just over lOOng . 
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Fig. 4. — Evolution of the mass-to-fiux ratio (in units of the critical value) fi c and of the Jeans number J c for the densest, moderately 
long-lived core in the subcritical run M10J4/3.01. The lower time-axis is given in units of the frame number in the animation. 
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Fig. 5. — Animations of the supercritical runs as in fig. 3, but with the density iso-surfaces corresponding to n = 10no (blue), n = lOOno 
(yellow), and n = lOOOno (red), a) Run M10J4/3.1. b) Run M10J4/31. The still frames in the paper edition show: (a) Frame 210 of run 
M10J4/3.1, at which time the third (and last) collapsing core (on the left) is forming. This core is fully collapsed by frame 220. (b) Frame 55 
of run M10J4/31, showing two of the long-lived but non-collapsing cores at the center (detailed in fig. 8) and left, and a collapsed core on the 
right. 




Fig. 6. — Iso-density surface 3D map of the clump-core system appearing at early times in run M10J4/3.1, detailed in fig. 7. The times 
shown are a) t = 0.06t s (frame 30 in the animation of fig. 5a ) and b) t = 0.08t s (frame 40). The iso-density surfaces are at n = 10no (blue), 
n = 40nrj (yellow) and n = lOOrao (red). The latter does not exist at frame 30. 
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Fig. 7. — Evolution of: (Left panel) the mass-to-fiux ratio (in units of the critical value) /j, c and of (right panel) the Jeans number J c for the 
first collapsing object in the moderately supercritical run M10J4/3.1. The various curves correspond to different values of the clump-defining 
threshold density m. Note that at frame 40 there are two cores defined by nt = 100no, described by the two dashed lines. At frame 30 the 
peak density is below 100no and therefore the dashed lines do not extend down to this frame. Note also that at frame 30 there are two cores 
defined by nt = 40rto, but their n c and J c are nearly idenstical, and so only one curve is shown. 
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Fig. 8. — Evolution of: (left panel) the mass-to-flux ratio (in units of the critical value) fi c and of (right panel) the Jeans number J c for the 
first long-lived rebounding clump in the strongly supercritical run M10J4/31. The curve coding is as in fig. 7. Note that this clump is always 
supercritical, but almost always sub-Jeans. The failure to collapse even though J c is slightly above unity for nt = 10no at frame 20 can be 
understood in terms of the non-accounting for turbulent support, and also slight imprccisions in our estimates of the core parameters (see 
text). 
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Fig. 9.— 
animation. '. 



Animation of the non-magnetic simulation M10J4/3oo. The still frame in the paper edition shows the final frame (#44) of the 
Note the larger number of cores in comparison to either supercritical run. 



